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ABSTRACT 

We derive the structural parameters of the recently discovered very low luminosity Milky Way 
satellites through a Maximum Likelihood algorithm applied to SDSS data. For each satellite, even 
when only a few tens of stars are available down to the SDSS flux limit, the algorithm yields robust 
estimates and errors for the centroid, position angle, ellipticity, exponential half-light radius and 
number of member stars (within the SDSS). This latter parameter is then used in conjunction with 
stellar population models of the satellites to derive their absolute magnitudes and stellar masses, 
accounting for 'color-magnitude diagram shot-noise'. Most parameters are in good agreement with 
previous determinations but we now properly account for parameter covariances. However, we find 
that faint satellites are somewhat more elliptical than initially found and ascribe that to the previous 
use of smoothed maps which can be dominated by the smoothing (roimd) kernel. As a result, the 
faintest half of the Milky Way dwarf galaxies {My ^ —7.5) is significantly (4(t) flatter (< e >= 
0.47 ± 0.03) than its brightest half (My < -7.5, < e >= 0.32 ± 0.02). From our best models, 
we also investigate whether the seemingly distorted shape of the satellites, often taken to be a sign 
of tidal distortion, can be quantified. We find that, except for tentative evidence of distortion in 
Canes Venatici I and Ursa Major II, these can be completely accounted for by Poisson scatter in 
the sparsely sampled systems. We consider three scenarios that could explain the rather elongated 
shape of faint satellites: rotation supported systems, stars following the shape of more triaxial dark 
matter subhalos, or elongation due to tidal interaction with the Milky Way. Although none of these is 
entirely satisfactory, the last one appears the least problematic, but obviously warrants much deeper 
observations to track evidence of such tidal interaction. 
Subject headings: Local Group — galaxies: dwarf 



L INTRODUCTION 

Owing to new sky surveys, the number of known 
dwarf galaxies or dwarf galaxy candidates that reside 
within the Local Group has d ouble d over the last 
very few years (IZuck er et al.l 120041: IWillman et al.l 



2005allbl: iBelokurov et al. 200 a iMartin et all 1200 



Sakamoto & Hasegawa 2006; Zuckcr et al. 2006alJ 



Belokurov et al 
2007; Walsh et al 
2OO81 : iMcConnachie et al 



2007; Ib ata et a l 
2007al: 



2007: Maicwski et all 



Zucker et al.ll2007l; llrwin et all 
I2008D . 



These new discoveries 
have extended the realm of galaxies to objects ^ 100 
times fainter than were known before. This enlarged 
population of extremely dim satellites is bringing us 
closer to an (at least) partial solution to the apparent 
overproduction of dark matter sub-halos in cosmological 
simula tions of Milky Way ( MW) and M31-like sys- 
tems (Koposov et ai." '2007a'; 'Peiiarrubi a et al.l l2008al : 
iSimon & Geha 2007; Strigari et al. 2003). This is 
also supported by the spectroscopic surveys of these 
objects and measurement of their velocity dispersions 
that require them to be among th e most dark matter 
dominated objects known to dat e jKlevna et al.' '2005'; 
Muiioz et al] 120061 : iMartin et all [2007; Simon & Geha 

Although some satellite can didates are still being 
searched for and proposed (e.g. iLiu et all |2008( ). most 
MW satellites within the reach of the Sloan Digi- 
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tal Sky Survey (SDSS ) have probably been discovered 
(jKoposov et al.l l2007af l . More and more effort is now 
spent on trying to understand whether these systems 
reflect the faint end of the dwarf galaxy luminosity 
function, whether they are disrupted versions of their 
brighter sibli ngs or whether they are entirely differ- 
ent objects (iBelokurov et all [200^ iGilmore et~aLll2007t 
iPefiarrubia et al.ll2008b| ). But as these studies are based 
on the observed properties of these systems, it is impor- 
tant to quantify their structural parameters or luminosi- 
ties rigorously, given their often low number of stars (a 
few tens to a few hundreds down to the SDSS flux limit). 

As a simple comparison, both Bootes I (Bool) and 
Draco (Dra) are systems containing mainly (if not only) 
metal-poor and old stars, at similar distances from us 
(~ 75 vs. 60kpc), yet their difference in luminosity 
{My = —8.8 vs. My — —6.0) means that the num- 
ber of stars observable in Bool is roughly 10% that of 
Draco under similar observing conditions. Neverthe- 
less, for the faint satellites, parameters are usually de- 
rived with the same techniques that are used for their 
brighter counterparts. In particular, the structural pa- 
rameters are mainly derived by assuming spherical sym- 
metry or through t he intensity-weighted second moments 
technique (see e.g. iMcConnachie fc IrwinI [20061 for a de- 
scription). But in faint galaxies, this technique requires 
some smoothing to produce maps that are not dominated 
by shot-noise and on which the method can be applied re- 
liably. The impact of any such smoothing scheme (pixel 
size, smoothing kernel size, threshold over which pixels 
are used) on the recovered parameters has not been in- 
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vestigated extensively. 

In addition, the 'total magnitude' of extremely faint 
stellar systems is conceptually not well defined since 
this parameter is usually determined from the luminosity 
of individual stars. This is perfectly valid when thou- 
sands of stars provide a good sampling of the color- 
magnitude diagram (CMD), but as the CMD becomes 
more sparsely populated in faint satellites, the recov- 
ered luminosity can be affected by the stellar evolu- 
tion state of individual stars, in particular how many 
of them are close to the tip of the red giant branch. 
In Willman 1 (Will) single stars could have a strong 
but transient impact on the luminosity of the system, 
currently estimated at ^ 1000 L0. Finally, the ex- 
tent to which the seemingly distorted shape s of most of 
these new syst ems [e.g. Boo tes I (Bool, Belokurov et al.| 
[2006[ ). CVnl (IZucker et all 2006b: Martin ct al. 2008), 
Coma Berenices (Com , Belokurov et al.n2007j ). Hercules 
(Her, EoEgian et al.l |2007D, Ursal^.ior II (UM all, 
IZucker et al1 l2006a^ and Will (IWillman et al.ll2006f )] are 
intrinsic or a reflection of the few stars available to map 
them also warrants close scrutiny. 

Therefore, we set out to re-derive the structural param- 
eters (centroid, half-light radius, ellipticity, position an- 
gle) and luminosities of the newly discovered MW satel- 
lites through a Maximum Likelihood (ML) algorithm. 
This algorithm has the advantage of directly fitting the 
best stellar density model to the data — the positions 
of individual stars — without any need for smoothing. 
Moreover, one of its outputs is the actual number of 
stars present in each system, down to the SDSS depth, 
which we can then use to determine its expected mag- 
nitude. This latter value then accounts for the scatter 
that is induced by the individual evolution of member 
stars. Finally, we use the best fit structural models 
to test whether shot-noise alone can explain the seem- 
ingly irregular shape of smoothed stellar density maps, 
or whether they deviate enough from these best models 
to require dumpiness in their underlying distribution. 

This paper is organized as follows: § 2 describes the 
Maximum Likelihood algorithm and presents its results 
applied to SDSS data. We also analyze discrepancies 
with previous estimates of the ellipticity of some sys- 
tems. In § 3 we study the impact of the sparseness of the 
systems on their absolute magnitude while the search for 
deviations from the best model is explained in § 4. Fi- 
nally, we use § 5 for a discussion of the ensuing homoge- 
neous dataset of morphological parameters and consider 
possible scenarios that could explain the higher charac- 
teristic ellipticity found in the faintest half of the Milky 
Way dwarf galaxies. § 6 concludes this paper. 

2. STRUCTURAL PROPERTIES 

2.1. Maximum Likelihood parameter estimates 

While the approach of estimating structural parame- 
ters for stellar systems directly from a star catalog via a 
Maximum Likelihood fit is in principle established (e.g. 
iKlevna et aDll99a iM^tfall et al.ll2006f ). it has not been 
used very widely. Therefore, we recapitulate the ap- 
proach here and spell out our particular implementa- 
tion. We start by presuming that the observables, in our 
case the positions of stars in the SDSS, have been drawn 
from a (spatial) model distribution that is described by 



j model parameters, pi,P2, ■ ■ ■ ,Pi- The g oal of the well- 
known ML technique (e.g. lLuptoniri993f ) is to find the 
set of parameters (pi,p2, ■ • ■ ,Pj") for which the observa- 
tions become most likely. This means maximizing the C 
function defined as: 



^(Pl,P2, ■ 



.Pi J 



1,P2, 



(1) 



where £i(pi,p2, . . . ,Pj) is the probability of finding the 
datum i given the set of parameters pi,p2, . . . ,Pj. The 
value of C is usually determined over a grid that explores 
the region of interest in the j-dimension parameter space 
in order to identify the global maximum C{pi,p2, . . . ,Pj), 
subject to a normalization constraint that has no im- 
pact in our analysis. In practice, £i(pi,p2, . . . ,Pj) is 
the value of the model defined by the set of parameters 
Pi,P2, ■ ■ ■ ,Pj for the properties of datum i. 

For the case of the MW satellites we are interested in, 
where the number of SDSS stars in each satellite can be 
very low (a few tens to a few hundreds over the studied 
region of tt deg^), we restrict the set of models for the 
member star distribution on the sky, Ss(''), to only ex- 
ponential profiles. They are specified by one parameter 
fewer than the more customary King profiles, are also 
easier to deal with analytically than the Plummer mod- 
els and usually yield reasonable fits^. We also account for 
the background level of stars, Sf,, assumed to be constant 
over the studied region of the sky. We therefore consider 
the model family E(r|pi,p2, . . . ,Pj) of stellar density in 
the SDSS, defined by the set of parameters pi,p2, . . . ,Pj, 
that can be expressed as: 



£i(pi,P2, . . . ,Pj) = S(r,|pi,p2, ■ ■ • ,Pj) 
= S'oexp(- — ) -I- Sfc 

re 

(2) 

with Sq the central stellar density, r,, the exponential 
scale radius of the profile and the elliptical radius co- 
ordinate of datum i, as we consider flattened models. We 
are now left with the task of relating these three param- 
eters (Sq, ri, re) to the five parameters that need to be 
constrained, namely: the centroid of the satellite defined 
by the two celestial coordinates ao and Sq or Xq and Yq , 
the spatial offset of the centroid compared to literature 
values; the position angle of the satellite from north to 
east 0; its ellipticity^ e; its size defined by the half-light 
radius^ . The background level Eh is the sixth parame- 

^ In fact, in most cases none of these profiles yields better results 
than the others as they do not account for the complexity that is 
being unveiled in an increasing number of MW satellites. Indeed 
several MW dwarf galaxies have now been shown to harbor multiple 
stellar populations (e.g. Sculptor, Fornax; see :McConnach ie et al.l 
120071 and references therein). However, they are a very useful and 
convenient way to parameterize and compare the structure of these 
galaxies. 

^ Throughout this paper, the ellipticity is defined ase = l — b/a 
with b the scale-length of the system along its minor axis and a 
that along its major axis. 

* The radius we determine in this manner is in fact the half- 
density radius but assuming there is no mass segregation in these 
systems, it also corresponds to the half-light radius that is com- 
monly used in similar studies. 
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ter constrained with the ML fit but its influence on S(r) 
is obvious in equation ([2|). 

A simple integration of Ss(r) out to the half-light ra- 
dius easily gives the first relation = 1.68re. Moreover, 
star i taken from the sample, with celestial coordinates 
(ofi. Si), has a spatial position {Xi, Yi) related to the cen- 
troid by^: 

( X^- Xq = (aj - ao) cos((5o) 

while the elliptical radius is related to the spatial position 
by: 



probes the region around the global maximum 
of C. The final grid has grid point sizes of 
(0.1', 0.1', 1.0°, 0.01, 0.1', 0.0001 stars/arcmin^). 

Finally, assuming the uncertainties on the various pa- 
rameters are well behaved near the global maximum (i.e. 
Gaussian or close to Gaussian), it can be shown that 
2 In (£(pi)) b ehaves as a distribution (e.g. § 10.3 of 
lLupton|[l993[ ) and the ka confidence intervals^ of the pa- 
rameters are bounded by the values that correspond to 
the function 2 In (£(pi)) dropping by k'^. In practice, this 
translates to relative likelihood values of, respectively, 
0.61, 0.136 and 0.011 for 1, 2 and 3a. 



n = (j^Y^^{X,cose-Y,sine)y+(^X,sme+YiCOse 

(4) 

With this expression, the major axis of the system is 
aligned with the position angle axis. 

Finally, 5*0 is constrained by the number of stars in the 
SDSS sample iVtot as this number must also correspond 
to the integration of the stellar density S over the region 
(of 1°) around the satellite, that is, with the integration 
performed on the circular radius r' from the centroid: 

iVtot= / / S(r')dr'r'd(/. 
Jo Jo 

= / S,(r')dr'r'd(/) + 36007rSb 

Jo Jo 
= N^ + 36007rSf,. 

(5) 

with N^: the number of stars in the sample that belong 
to the satellite and S, and St, expressed in units of 
stars/arcmin^. Given that, for reasonable values of r^, 
T,s is negligible at r' = 60', one can integrate the first 
term to r' — +oo instead. From there a change of vari- 
ables from the circular radius r' to the elliptical radius r 
used to express in equations ((2]) and fl]) is required to 
simplify the integration. Its Jacobian is (1 — e), yielding: 



271- 



^0 



S,(r')dr'r'. 



2t! 



(1 - e) / / Y.s{r)drrd<j) 
Jo 

2 



= (l-e)2^5or 
from which we eventually derive 

JO 



1/2 2.2. Results 

We then proceed to apply the algorithm to consis- 
tently derive the structural parameters of the 12 Milky 
Way satellites recently discovered in the SDSS [Bool, 
Bootes II (BooII), CVnl, Canes Venatici II (CVnII), 
Com, Her, LcoIV, LeoT, Segue 1 (Segl), Ursa Major I 
(UM al) , UMall and Will], the 5 candidates proposed by 
Liu et al.l (|2008[1 and Objects X, Y and Z proposed by 



(6) 



(7) 



2^r2(l-e)- 

By substituting r^, Vi and Sq in equations (H]) and ^ 
by their expressions, C is now directly expressed as a 
function of (ag, i5oj 6', e, r/j, E^). 

The set of best parameters (ao, ^O: e, ''I, ^h) is 
determined through an iteratively refined grid that 

^ It should be noted that this expression is only valid for stars 
relatively close to oq and (5o, as is the case in this study. Large 
deviations from the centroid would require properly accounting for 
the projection of stars on the tangential plane of the sky. 



Koposov et al.l ()2007a^ . To test the quality of its results. 



we also apply the algorithm to the well studied Dra dwarf 
galaxy 

To do so, we extract stellar objects from th e SDSS 
Data Release 6 (jAdelman-McCarthv et al.ll2008( ) m a re- 
gion of radius 1° around the center of the satellites. We 
then apply a selection box loosely constructed by eye 
around their CMD features ((sub-)giant branch, horizon- 
tal branch, main sequence turn-off and main sequence in 
the closest cases) so as to exclude clear Milky Way fore- 
ground or background contaminants. Additional mag- 
nitude cuts at r < 22.0 and g < 22.5 ensure that the 
photometry is of good quality. I n the case of the very 
faint objects Will as well as the iLiu et all (|2008D and 
iKoposov et al.l ()2007a[ ) candidates, we push these lim- 
its to r < 22.5 and g < 23.0 in order to increase the 
number of stars. For each satellite, the resulting stellar 
spatial distribution is then used to determine the best 
stellar density exponential model through the ML algo- 
rithm described above. 

The algorithm converges, i.e. finds a unique likelihood 
maximum, in all cases ex cept for 4 of thelLiu et all (j2008l ) 
candidates and 2 of the lKoposov et all ()2007aD ~obiects. 
despite the deeper magnitude limits (see below for more 
details). The results of the fits are summarized in Table 1 
where the best values of the six parameters are listed 
(ao, Sq, 6, e, and A^*, the integrated number of stars 
in the galaxy down to the assumed SDSS depth). The 
resulting profiles are built in Figure [U and compared 
to the binned data assuming the best elliptical model 
that is found. In all cases, the two-dimensional relative 
likelihood distributions are shown in the Appendix with 
the contours representing drops of 50, 90 and 99% in 
likelihood compared to the best model. Most parameter 
combinations exhibit no significant covariance, except for 
r/i, e and Eb. As expected, satellites with more stars 
yield better constrained models (e.g. Bool, CVnl, Dra) 
but even sparsely populated systems show well behaved 

^ The confidence interval is here defined in the usual sense of 
a normal distribution as the regions that contain the central 68.3, 
95.4, 99.73% of the probability distribution function in the case of 
1, 2, 3(T respectively. 
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likelihood distributions. The exception is the position 
angle 9 in some cases, but this is expected as for (quasi- 
)round systems, this parameter is ill-defined (e.g. BooII, 
LeoIV). 

Below follows an object by object discussion of the 
derived parameters. 

Dra - Applying our ML fitting to the Draco galaxy 
provides a good opportunity to verify that the algorithm 
works properly on a galaxy that has been well studied. 
Our structural param eters are very s imilar to previous 
determinations bv Odenkirchen et al.l (|2001al 6 = 88 ±3° 
and e = 0.29±0.02), also from SDSS data. The half-light 
radius we measure is in excell ent agreeme n t with that 
obtained from deeper data by ISegall et al.l (|2007l . = 
10.09 ± 0.02', although this value is derived from their 
best Plummer model) . This comparison gives confidence 
in the results of the ML algorithm. 

Bool - All the parameters derived for this galaxy 
are in agreeme n t wit h previous measurements from 
iBelokurov eFall (|2006l . 61 = 10 ± 10°, e ~ 0.33 and 
rh = 12.6 ± 0.7'). We derive a slightly higher elliptic- 
ity but, given the uncertainties, our value is comparable 
to e ~ 0.33 (but see below for a detailed comparison). 

BooII - No determination of structural parameters 
exists for this puny object except for a h alf-light ra- 
dius o f 4.0 ± 1.9' in the discovery paper (Wals h et al.l 
l2007al) and 3 . 6 ib 1.2 ' from deeper follow-up observations 
(| Walsh et al]|2007b[ ). The ML fit gives the picture of a 
rather spherical system. 

CVnl - As for Bool, this galaxy is well populated in 
the SDSS, and the ML values are all within la of those 
of the discovery paper (jZucker et al.ll2006bl = 73 ± 3°, 
e ~ 0.38 and = 8.4 ±0.5'). 

CVnII - This is a faint and distant object that does 
not have many stars within the SDSS, which proba- 
bly explains why our parame ters deviate from those of 
Sakamoto &: Haseeawal ()2006L e ~ 0.1 and — 0.9') or 
Belokurov et all (|2007l . 6* ~ 0°, e ~ 0.3 and rh 3.3'). 



However, since none of these authors quote any uncer- 
tainties, it is difficult to ascertain if these differences 
could be due to poorly constrained values. The rela- 
tive likelihood contours for this satellite (Fig. [U 4*'' set 
of panels) nevertheless exhibit gentle variations and we 
have no reason to expect the ML algorithm to have failed 
and be stuck on a local maximum of the parameter space. 
Our best model also yields a stellar density profile that 
is very close to the data (Fig. 1). We find that CVnII is 
a rather elliptical system with an axis ratio of ~ 2 : 1. 

Com - With this reasonab ly populated obj e ct, we are 
once again in agreement with lBelokurov et al.l ()2007L 9 ~ 
—60°, e 0.5 and rh — 5.9'). We may determine the 
system to be slightly rounder, but e ~ 0.5 is just within 
our uncertainties. 

Her - This galaxy is one of the recently discov- 
ered galaxies that deeper follow-up data from the Large 
Bin ocular Telescope (LB T) revealed to be more ellipti- 
cal (jColeman et al.ll2007 . 9 ~ -7 3°, e - 0.65) than from 
the smoothed maps of the SDSS (jBelokurov et al.l [20071 
9 ~ —60°, e ~ 0.5). Since it could be the sign that 
this latter value is driven by the smoothing of the shal- 
low data, it is particularly interest ing that our ML value 
(e = 0.68t°;o8) is very close to the C oleman et al.l ()2007( ) 
value. Even from the SDSS data alone, we confirm that 
Her is one of the most elliptical satellites of the MW. On 



the other ha nd, the half-light radius we measure is simi- 
lar to that of IBelokurov et al.l ()2007L rh ^ 8.4'). The one 
measured from the LBT dataset was erroneously listed as 
4.4 ± 0.3' where it should have been 9.4 ± 1.4' (Matthew 
Coleman, private communication), once again very sim- 
ilar to the value of the ML fit. 

LeoIV - This almost spherical object contains very few 
stars in the SDSS (A^* = 30 stars), leading to sizeable 
uncertainties. We are therefore in good agreement with 
IBelokurov etlll (pOOTl . 9 ~ -5°, e ~ 0.25 and ~ 3.4') 
although our fit converges towards a somewhat smaller 
system. 

LeoT - Sitting in the outskirts of the MW at ~ 
400 kpc, this peculiar galaxy has unti l now managed to 
retain some gas in the form of Hi (| Irwin "erain2007l : 
iRvan-Weber et all [2008,1 and has undergone a signif- 
ica nt episode of sta r formation over the last Gyr or 
so (|Irwin et al.l[2007t Ide Jong et al.ll2008aD . The SDSS 
data exhibit both young and old stars that are re- 
vealed to_Jonow_differeiit density profiles from deeper 
data (jde Jong et all l2008af) but given the low number 
of LeoT stars in the SDSS, it is necessary to merge both 
young and old stellar populations in the same sample 
to allow the ML algorithm to converge. The parame- 
ters listed in Table 1 are therefore the overall parameters 
of the dwarf. They are neverthe l ess rea sonably close to 
those determined by llrwin et al.l (I2007L 9 :^ 0°, e ~ 0.0 
and rh ~ 1.5') from SDSS data or lde Jong et al.l (|2008al 
rh — 1.0 ± 0.1') from deeper LBT data. In this latter 
case, the difference in the rh value could also be due to 
the different depths in the two datasets and mean that 
the shallow data are dominated by the more extended 
old population (r^ = 1.2 ± 0.1' from the LBT data). 

Segl - Despite being one of the faintest MW satel- 
lites, the proximity of Segl to the Sun ensures it contains 
enough stars in the SDSS for a good fit. As for Her, we 
find the system to be more elliptical and slightly rotated 
when c ompared to the values determined from smoothed 
maps (jBelokurov et all l2007l . 6* ~ 60°, e ~ 0.3). The 
half-light radius is however in excellent agreement with 
their value (r^ ~ 4.6'). 

UMal ~ Deep data down to the main-sequence turn- 
off of this galaxy obtained bv lOkamoto et al.1 (|2008f l give 
structural parameters {9 ~ 78° and rh — 11.3') that 
are very similar to those yielded by the ML algorithm. 
We confirm the larg er size of the system compared to 
the discovery paper (jWillman et al.ll2"005E rh 7.75') 
but , once again, derive a higher ellipticity than the 
one lOkamoto et all (|2008D estimate from their smoothed 
maps (e ~ 0.54). Our value makes UMal the most flat- 
tened MW satellite galaxy. 

UMall - This is yet another case were we derive 
an ellipticity that is high er than that determ ined from 
smoothed maps (e ~ 0.5; iZucker et al.l[2006ah . We are 
also able to refine the value of rh that is within the large 
range initially given by these authors. 

Will - Table 1 gives the first measurements of the el- 
lipticity and position angle for this satellite, once more 
revealing a rather elliptical system although not as dras- 
tically so as UMal, UMall or Her. The ML best fit 
has a half-light radius similar to p revious estimates from 
deeper d ata (r/, - 1.75 ± 0.5', IWillman et al] l2005al : 
rh ~ 1.8'. [Martin et al.ll2007f) . 

SDSSJ08U+5105, SDSSJ0821+5608 & 
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TABLE 1 

Derived properties of the satellites 





Bool 


BooII 


CViiI 


CVnII 


Com 


Qo (J2000) 


+ 14''00"03.6'' ± 1.5" 


+ 13''58"08.0'' ± 2.0" 


+ 13''28"03.9° + 1.3° 


+ 12''57"10.0° ± 0.4° 


+ 12''26"57.8° + 2.5° 


i5o (J2000) 


14°30'42" ± 43" 


12°50'54" ± 34" 


33°33'33" + 13" 


34°19'33" ± 13" 


23°55'09" + 19" 


e (dog.) 


+ 14 ± 6 


-35+f, 

— 55 


+70^4 


-3+9 


-65+^" 


(arcmin) 


n 3Q + n Dfi 
12 6+^" 


0.21 ± 0.21 

^■■^-1.4 


39 + n 03 
8.9 + 0.4 


52+'''^° 

^^•"-0.2 


"■'"'-0.14 

6.0 + 0.6 


N, 


3241^1 




463 + 18 


291= 


99 + 13 


E(B - V) 


0.017 


0.031 


0.013 


0.010 


0.018 


D(kpc) 


66 ± 3'"' 


42 ± 8<''' 


218 + lO'"' 


160+*(<=' 

— 5 


44 + 4<'^' 


1^1. ( T» 1 
' h \ P<-J 




51 ± 17 


564 + 36 


'*-10 


77 + 10 


My 


fi -t- n 9 


9 7 + n Q 


^8.b_(, 1 


A 4- n ^ 
— IE u.o 


/I 1 + n ^ 


^scatter (m^'g-) 


±0.2 


±0.8 


+0.1 


±0.5 


+ 0.4 

-0.5 






28 1 ± 1 6 


27 1 + 2 


26 l+°-'' 


97 O + 0.7 


Lv(Lo) 


3.0 ± 0.6 X 10* 


1.0 ± 0.8 X 10^ 


2.3 + 0.3 X 10^ 


7.9I3 ;* X 10^ 


3.7i;;6 X 10^ 


M.,Kroupa(M0) 


3.4 ± 0.3 X 10* 


i.4j;!;-5 X 10^ 


3.0 + 0.2 X 10^ 


8.0lj-° X 10^ 


4.8 + 0.9 X 10^ 


M..Salpotcr(M0) 


6.7 ± 0.6 X 10* 


2.81^ X 


5.8 + 0.4 X 10^ 
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Note. — All struct u ral pa rameters {ao, So, 6, e, rf^ and A^*) are determined through the M L fit. The extinction value s are derived 
from th e ISchlegel et"all ((T998I) maps and the distances are taken from the following papers: (a) IDairOra et al.l f200^, (b) IWalsh et al] 
(I20q7bl). (c)IMartiri et al. (2008), (d) Bclokurov ct al'.' (2007), (c) Greco ct al. (2008), (f) Coleman ct al. (2007), (g) dc Jong ct al. (2008a!)7 
(h) lOkamoto et al.l ( 12008 ). (i) Zucker et a l. (2006a), (j) Willman ct al.. (2005a ). (k) Liu ct al. (2008i) . (1) ,Bonanos ct al.. (2004 ). The 
magnitudes of the systems, My, their scatter induced by 'CMD shot-noise', (Jscatter, their central surface brightness, fio v, ^^d their stellar 
masses, M^ Kroupa and Salpeton a-i's ^ll determined as explained in § 3. 

(t) We refrain from determining the magnitude of LeoT given the presence of multiple stellar populations that render its determination 
very uncertain (see § 3 for more details). 

(I) In the absence of any distance estimate for Object Y we cannot derive parameters other than those directly estimated through the ML 
fit. 
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Fig. 1. — Stellar profiles of the dwarf galaxies. In all panels, the stellar density measured in fixed elliptical annuli from the SDSS data 
is shown as dots, using the best ellipticity, position angle and centroid found by the ML algorithm (the uncertainties are derived assuming 
Poisson statistics). The full lines represent the best models constructed from the parameters of Table 1. They are not fits to the data 
points but show very good agreement with them, even in poorly sampled satellites such as e.g. BooII, LeoT or Segl that all contain less 
than 100 stars. 
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SDSSJ1329+2841 - The ML fits do not yield sig- 
nificant resul ts for thes e three MW satellite candidates 
proposed by iLiu et all (|2008l ). even when using mag- 
nitudes cuts 0.5 mag fainter. They correspond to 
overdensities of only 3.5, 2.0 and 2.5ct in N^, respectively 
and for the last 2 objects the best ellipticity found by 
the algorithm is an unphysical e = 1.0. Even though 
we cannot rule out that the ML is just not sensitive 
enough, none of these objects is visible on maps of the 
spatial distribution of SDSS stars, contrary to all the 
other objects (even the extremely faint Segl, Will or 
BooII). 

SDSSJlOOO+5730 - From the inspection of the SDSS 
catalog around this object, it appears that this "overden- 
sity" is in fact within ~ 25' of a series of holes in the pho- 
tometric catalog that result frqni bad seeing. These holes 
seem to have driven iLiu et all ([2008) to tag this under- 
density of the background (that is an order of magnitude 
smaller than for their other candidates in their Table 2) 
as an over density produced by a MW satellite. 

SDSS J 1058+2 843 - The last candidate of the 
iLiu et all (|2008[ ) list is the only one for which our al- 
gorithm converges, revealing a very sparse object with 
sizable uncertainties on its parameters. Our measure of 
the half-light radius is 1.5cr smaller than in the discov- 
ery paper but this is not surprising given the scantness 
of this system. Deeper photometry and/or spectroscopy 
are required to definitely confirm that this overdensity 
is indeed a MW satellite. But since its properties are 
similar to those of Segl and Will, we will still consider 
it as a satellite in our analysis (removing or keeping this 
satellite does not impact our overall conclusions). 

Objects X & Z - As for the three first candidates of 
iLiu et all (|2008f l. the algorithm does not converge on 
physical va lues for these two candidates proposed by 
iKoposov et al. (2007al). 

Object Y - This is the one object of the three 
iKoposov et al.l (|2007af ) candidates for which the algo- 
rithm converges. Although the derived values are quite 
uncertain, the data and best profile agree reasonably well 
in Figure 1 and give the picture of a rather elliptical 
ove rdensity. This object is unfortunately too sparse for 
the lde Jong et al.l ()2008b| ) CMD fitting technique to con- 
verge and the ensuing absence of any distance estimate 
prevents us from deriving the physical size of this object. 

In summary, the ML algorithm behaves very well as is 
confirmed by the comparison between the preferred pro- 
files and the binned data in Figure [H Although there 
are a few cases where some scatter is present, there is 
usually a very good correspondence between the radial 
profiles. Where available, the derived structural param- 
eters are also in good agreement with previous measure- 
ments. It is particularly reassuring to notice that in the 
case of the well populated D raco galaxy, our pa r ameter s 
are very similar to those of lOdenkirchen et al.l (|2001af ). 
In the few cases where deeper follow-up was obtained on 
the less populated faint satellites, the ML algorithm ap- 
plied to the shallower SDSS data yields results that are 
consistent within the uncertainties. This is the case for 
the ex tremely faint BooII (comparison with lWalsh et al.) 
I2007hl ^ or the very elongated He r, for which we derive 
again an axis ratio of 3:1, as in IColeman et al.l ()2007l ) 
from LBT data reaching ^ 4 magnitudes deeper. 

The case of Her is also a good example to illustrate how 



smoothing the maps before deriving parameters through 
the intensity weighted second moments method appears 
to lead to measurements of smaller ellipticities: ~ 0.5 
vs. 0.68tg |]^ for Her but also ~ 0.54 vs. 0.80 ± 0.04 for 

UMal, ~ 0.3 vs. 0.521°;}° for CVnH, ~ 0.3 vs. 0.48+^;}!] 

for Segl or ~ 0.5 vs. 0.63l^;[55 for UMall. 

To investigate this effect closer, we build 500 mock 
models of galaxies that are similar to Bool, one of the 
typical new satellites. They all have r/j = 12.6', iV* = 
324 stars and an ellipticity, einput that we fix. To re- 
main as close as possible to the observed SDSS data, we 
also add randomly scattered background with the den- 
sity found by the ML algorithm. The colors and mag- 
nitudes of these contaminants, that are chosen to fall 
within the selection box used to derive the properties 
of Bool, are randomly picked from all stars in the SDSS 
within a region of 5° x 5° around Bool, excluding the cen- 
tral 2° X 2° region. We then recover the ellipticity of the 
models, eoutput, using both the intensity- weighted second 
moments method and the ML algorithm. In the case of 
the former, a map of each system is built with pixels of 
2' X 2' that is smoothed with a round kernel of 4' disper- 
sion. Only pixels that are at least 3cr over the mean back- 
ground level are used (the mean background level and its 
r.m.s., fj, are determined from the circular annulus cen- 
tered on the model and with 0.75° < r < 0.9°). For 
the latter, we automate the ML algorithm that searches 
through an iteratively refined grid. We have checked on 
the Bool SDSS data that the results obtained through 
the iteratively refined grid are very similar to those listed 
in Table 1 for that object. The recovered ellipticities, 
Coutput, are compared to einput in Figure [2] for both meth- 
ods. The usefulness of the ML technique is manifest as 
it recovers the input ellipticity in almost all cases. Us- 
ing smoothed maps, on the other hand, leads to system- 
atically lower recovered ellipticities and this effect gets 
worse as einput increases. Both methods have issues for 
the rounder systems whose recovered ellipticity is slightly 
higher than the input one; this effect is due to the degree 
of freedom afforded by 9 for the (quasi-)round systems 
in the presence of a few outer stars. 

This example shows the limitations of smoothed maps 
in determining the structural parameters of faint satel- 
lites. Not only do they require some fine-tuning to have 
a set of pixel size, kernel size and background threshold 
that translate into a suitable map, but even with such a 
map, the ellipticity is measured rounder than it truly is. 
This effect becomes even more important as the number 
of stars in the system shrinks. The strongest deviations 
are indeed observed for systems that are less populated 
than Bool: CVnII, Her, UMal or Segl. Since in most 
cases, this e is then used to construct the radial profile 
of the satellites, this can have a significant impact on 
the derived structure of the system. In contrast, the ML 
fit is directly applied to the data (the spatial distribu- 
tion of stars), and all the parameters are derived at the 
same time, leading to best values and uncertainties that 
include the infiuence of the determination of the other 
structural parameters. 

While this analysis of the behavior of the ML fit is only 
done in detail for Bool-like systems, it should be noted 
that the efficiency of the algorithm actually scales with 
the density of stars in the satellite and not solely with 
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^input 

Fig. 2. — Comparison of the input ellipticity of the models, einputi 
with the output ellipticity, Coutput , recovered through the intensity- 
weighted second moments method applied to smoothed maps (filled 
circles) and through the ML algorithm (hollow circles) for Bool-like 
systems. The points represent the median value of 500 models and 
the error bars on the recovered values represent the central 68.3% 
of the distribution. To guide the eye, the full line represents the 1- 
to-1 correspondence that should be followed if the methods have no 
systematics: this is clearly not the case from the smoothed maps as 
toutput is systematically lower than einput. The dashed-dotted and 
dashed lines show the quality of the ML fit for lower density UMaL 
and UMall-like systems respectively. Uncertainties are similar to 
those of the Bool-like models but are not shown to avoid cluttering 
the Figure. In both cases, the recovered ellipticity falls in the region 
where the algorithm follows the 1-to-l correspondence (eoutput = 
0.80 and eoutput = 0.63 respectively). 

their number of stars N^, within the SDSS. Since a lower 
hmit of the density p can be estimated by assuming a 
round morphology, it scales as p oc ^ . A look at Table 1 

shows that p is higher than that of Bool in all other 
cases, except for UMal, UMall and Object Y. Therefore, 
apart from these three systems, the quality of the ML fit 
remains at least as good as shown by the hollow circles 
in Figure [2] and only very round systems can be found to 
be slightly more elliptical than they really are. The cases 
of the less dense UMal and UMall are also shown in the 
Figure as the dashed-dotted and dashed line respectively 
(the dashed line also applies to Object Y as its density is 
similar to that of UMall). Even though the overestimate 
of small ellipticities is more significant than previously, 
the ellipticities of these three systems are so high that 
they should not be hampered by the algorithm: in both 
cases, the recovered ellipticity falls in the region where 
the algorithm behaves properly and follows the 1-to-l 
line. 

Generally, the over-prediction of the ellipticity of very 
round systems could become a difficulty in the case of 
systems sparser than those considered here. But re- 
cently disc overed satellites all lie at the detection limit of 
the SDSS (|Koposov et al.|[2b07af) and most fainter (and 



sparser) candidates, s uch as the iLiu et all (|2008f ) and 
iKoposov et al.l (|2007ar ) ones, are not dense enough to al- 
low the algorithm to converge. Their confirmation def- 
initely requires deeper datasets than the SDSS, leading 
to denser samples on which to apply the ML technique. 
This will consequently elevate the issue of the overesti- 
mation of small ellipticities. 

3. LUMINOSITIES AND STELLAR MASSES 

The use of the total magnitude to quantify the amount 
of stars in these very sparsely populated systems is prob- 
lematic. Their observed luminosity can indeed vary 
strongly through time, depending on whether a few stars 
have evolved along the red giant branch or most of 
them remain along the main sequence or main sequence 
turnoff. To remove the influence of this 'CMD shot- 
noise', we propose to characterize the systems by their 
number of member stars (above a luminosity threshold) , 
N.^,, instead of the luminosity of these stars, and study 
the impact of N,, on the total magnitudes and stellar 
masses that are likely to be observed. 

For each satellite^ we have stellar population m odels, 
also derived from SDSS data (jde Jong et al.l[20"08bf ) . and 
the ML fit gives A^* with its associated uncertainty; so 
we can generate numerous CMDs of each system. The 
luminosity of the stars in each of these CMDs then 
yields the 'observed magnitude' of this particular real- 
ization of the satellite. We then define the 'expected 
absolute magnitude' of the system as the median of 
the distribution of observed magnitudes®. In detail, a 
high signal-to-noise CMD is created for each dwarf, con- 
taining on average 70 000 stars, which includes a realis- 
tic ph otometric error and c ompleteness model for SDSS 
data ()de Jong et al.ll2008bf l. These 'master-CMDs' are 
populated assuming a standard Salp eter or K roupa ini- 
tial mass function (Salpeter 1955; Kroupa et al. 199^ 
down to the hydrogen burning limit (0.08 Mq) and a 
binary fraction of 0.5^. Since all the satellites consid- 
ered here are dominated by old populations, the stellar 
ages are spread between 10 and 13 Gyr in all cases and 
the metalliciti es are chosen to corre spond to the best- 
fit values from Ide Jong et all (|2008bl . given in their Ta- 
ble 5), where available. As these metallicities are only a 
means to parametrize and reproduce the SDSS CMDs, 
we abstain from using spectroscopic metallicities. De- 
spite being more accurate, these could be systematically 
offset from the isochrones and lead to artificial CMDs 
that are different from the observed ones. In the case 
of BooII, Draco and SDSSJ1058-h2843, whose stellar 
population models were not derived by Ide Jong et al.l 
(j2Q08b), we assume the same old populations and metal- 
licities of [Fe/H] = -2.0 (I Walsh et al.l l2007;j) . -1.8 
(jAparicio et al.ll200lD and -2.0 ()Liu et'aLll2008D respec- 
tively. In order to track not only the luminosities of the 
satellites, but also their stellar masses, the mass of each 

^ Except for Object Y, whose lack of a distance estimate prevents 
us from determining its absolute magnitude. 

* As these distributions are close to being normal distributions, 
this is almost equivalent to taking the mean of the distributions. 

^ In Draco, IKlevna et al.l II2002I ) have shown that there is mild 
evidence of a lower binary fraction than the solar neighborhood 
value used here. However, assuming a fraction of 0.3 does not 
significantly change our derived magnitudes: in the case of UMall 
which, owing to its proximity, has one of the deepest CMDs of all 
the satellites, the measured magnitude changes by less than 0.05. 
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star in the 'master-CMDs' is stored, together with its 
magnitude and color. 

One thousand realizations of each satellite are gener- 
ated randomly by selecting TV, stars (with its associated 
uncertainty listed in Table 1) from the 'master-CMD' 
to fall within the selection box previously used to select 
member stars. The 'observed magnitude' of each realiza- 
tion is then determined from the luminosity of all stars 
above a limiting magnitude of g = 21.5 (half a magnitude 
brighter than the approximate 100% completeness limit 
of the SDSS), transformed to the F-band using color 
equations from iJester et al.l ()2005[ ). Using the distance 
listed in Table 1, the luminosity of stars below the magni- 
tude limit is accounted for by using luminosity function 
LF) corrections derived from the LFs of iDotter et al.l 
2007f ) with the metallicity and age closest to each model 
population. 

In parallel, we also determine the stellar mass of the 
satellites by determining the mass down to the same mag- 
nitude limit as above and then correcting for the 'unob- 
served' stars present in the 'master-CMDs', down to the 
hydrogen burning limit. This correction^'^ is done for 
both a Salpeter and a Kroupa IMF to yield two esti- 
mates of the stellar mass listed in Table 1 for each satel- 
lite (M,,saipeter and Af,, Kroupa respectively). 

The resulting distributions of these 'observed magni- 
tudes' and 'observed stellar masses' are shown in Figure[3] 
along with their 'expected absolute magnitude' and 'ex- 
pected stellar mass', taken to be the median value of the 
distributions. The scatter produced by the 'CMD shot- 
noise' is readily visible from the horizontal error bars that 
represent the central 68.3% of the distributions. All val- 
ues and their zblcrgcattor are listed in Table 1. As could 
be expected, bright galaxies such as CVnl and Dra have 
an extremely small scatter, as their iV* is very well deter- 
mined by the ML fit. N^, is also large enough that there 
is a good sampling of the CMD: whatever the random 
realization, the resulting CMD is very similar. On the 
other hand, as the satellites become fainter and fainter, 
their CMDs are only sparsely sampled for each realiza- 
tion and the resulting magnitude is strongly influenced 
by the presence of only a few very bright stars. The stel- 
lar masses are much less affected by these CMD sampling 
problems; indeed the evolution of an individual star does 
not change its mass as strongly as its luminosity. 

For a proper comparison of our absolute magnitudes 
with the literature, we need to include the uncertain- 
ties that come from our assumptions on the age and 
metallicity of the dwarf and thence on the LF correc- 
tion. This correction is, however, almost insensitive to 
±0.5 dex changes in the metallicity (less than 1% change) 
while a shift of ±2 Gyr yields small changes of only ~ 3% 
in all cases. This translates into uncertainties of less than 
clf = 0.06 mag. Another source of uncertainty can come 
from the assumption that the satellites are well repre- 
sented by a single population model. In fact we know in 
the case of CVnl that this is not the case as ^ 5% of the 
stars are found to be young (jMartin et al.ll2008h . How- 
ever, such a low proportion of young stars only minutely 
affects the luminosity correction (by less than 1%) and 

In the stellar mass range probed by the SDSS photometry of 
these old populations, the slope of both IMFs is the same, so that 
this extrapolation is valid in both cases. 



does not impact our results. This is not the case anymore 
for systems with a significantly extended star formation 
history and we therefore refrain from determining a mag- 
nitude for LeoT. 

Finally, we also account for the uncertainty coming 
from the distance measurements, Udist (usually of the or- 
der of 10% in distance or 0.3 mag; see Table 1). The 
final uncertainties given in Table 1 correspond to the 
quadratic sum of all three uncertainties: (Cgcattcr + "'lf ~^ 
(T^jgj)^/^. Comparison of these with values of Cscatter 
readily shows that the 'CMD shot-noise' is the main 
driver of uncertainties on the 'expected absolute mag- 
nitude' of such faint systems. 

In the end, almost all of our luminosity measurements 
are less than la from literature values but with slightly 
smaller error bars. The D r aco v alue is in excellent agree- 
ment with that of lMateol (|1998l . My = -8.8 ± 0.2). We 
confirm the extremely low luminosity of BooII, UMall, 
Will and Segl, although for this latter case, we find it 
is even fainter than previously measured (— 1.5^Q g vs. 
—3.0 ± 0.6). At the other end of the spectrum, we find 
that CVnl is somewhat brighter than previously inferred, 
which makes it almost as bright as Draco: this galaxy 
should probably not be seen as a peculiar faint system 
but mainly owes its late discovery to its large distance. 
Our value for UMal is significantly smalle r than that of 
the discovery paper (jWillman et al.ll2005U My ^ -6.75) 
but, as for its^half-light radius, it confirms the new mea- 
surement by Zuckcr (200^ My — —5.5). This value also 
solves the issue of Bool containing more sta rs than UMal 
while being an apparently fainter system (jMufioz et al.l 
[20061) : Bool is brighter than UMal and should contain 
more stars. 

From the luminosity and structural parameters of the 
systems, we also estimate their central surface brightness 
Hoy, listed in Table 1. As the central surface bright- 
ness (expressed in units of LQ/arcsec^) is proportional 
to (1 — e)~^ [see equation (6)], our newly derived pa- 
rameters sometimes have a strong impact on the derived 
fJ-o.v that can be significantly different from previous es- 
timates. When our morphological model is similar to 
previous models, we nevertheless reach very similar val- 
ues (e.g. Dra with fioy = 25.0 ± 0.2 mag/arcsec^ vs. 
Hoy = 25.3 ± 0.5 mag/arcsec^ in iMateol fl99l or Bool 
with Hoy = 27.5 zfc 0.3mag/arcsec^ vs. Hoy = 27.8 ± 
0.5mag/arcsec^ inlBelokuro v et ai]|2006[ ). We find that 
there is no system fainter than Hoy — 28.0 mag/arcsec^, 
a threshold that is brighter than previously claimed but 
in very good agreeme nt with the detection limit of the 
iKoposov et al.l (j2007a[ ) automatic search. 

Apart from the scatter due to the uncertainty in N^, 
the largest contribution to the uncertainty of any stel- 
lar mass determination comes from the assumed IMF. 
The strength of this effect is readily visible by compar- 
ing the two stellar mass estimates listed in Table 1 for 
the Salpeter and Kroupa IMFs. 

It should be noted that even in the case of very faint 
satellites, the changes in the luminosities of the sys- 
tems have only a small impact on th e determinations of 
their total (dark an d luminous) ma ss (iKleyna et al.ll2005t 
iMuiioz et al.' '20061 iMartin et all 120071 : ISimon fc Gehal 
2007; Strigari ct all 120071 ^ as these are very model- 
dependent and strongly dominated by other uncertain- 
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Fig. 3. — Impact of the 'CMD shot-noise' on the 'observed magnitude' and 'observed stellar mass' of the satellites. Each histogram 
corresponds to a different galaxy and shows the distributions of 'observed magnitudes' (left) and 'observed stellar mass' (right) for different 
realizations of their CMD. The median value of each distribution, the 'expected absolute magnitude' (resp. 'expected stellar mass') that 
best characterizes the luminosity (resp. mass) of the system, is shown as a filled dot and the horizontal error bar represents their lo-gcattor 
(the region corresponding to the central 68.3% of the distributions). As expected, faint satellites show a much larger scatter in their 
magnitude since the presence/ absence of a few very bright red giant branch stars will strongly impact the recovered value. The masses 
show a smaller scatter as they are less sensitive to the evolution of the stars. 



ties. 



4. ARE THE STELLAR DISTRIBUTIONS 
'DISTORTED'? 

The best-fit model of a satellite can also be used 
to test whether the data imply deviations from this 
best (smooth and symmetrical) model. One of the 
striking features of most of the recently discovered 
systems is their apparent dumpiness, even when the 
stellar distribution is smoothed. It has been sug- 
gested on numerous occasions that their distorted 
isophotes could be the sign of tidal dis t ortion or in- 
trinsic dumpiness (iBelokurov et all 120061: IZucker et al.l 
2006at IW illman et al.' '2006': 'Bclokurov et all l2007t 
Coleman e t al. 2007; Martin ct al. 2008). However, the 
(very) low number of detected stars may leave the same 
impression, even if the parent distribution were smooth 
and symmetric. A good visual example of such a 'low 
number statistic' behavior is given in Figure 21 where 
the smoothed map of Bool (left) is compared to simi- 
lar maps built from purely spheroidal models with the 
structural parameters obtained from the ML algorithm 
e = 0.39 ± 0.06, = 12^6^^:^ and = 3241^8 (the 
background contamination is accounted for in the same 
way as in the models presented in § 2). It is readily ap- 
parent that, even with the smoothing that tends to drown 
small clumps, the realizations of purely spheroidal mod- 
els appear just as distorted as the observed data. Taken 
at face value, most of them would in fac t be sus pected 
to contain 'tidal features'. I Walsh et al.l ()2007bf) reach 
similar conclusions in their analysis of BooII. 

To better estimate the (non-) significance of the appar- 
ent distortions, we determine, for all satellites, the frac- 
tional r.m.s. deviation dsc/total of the data compared to 
the best model, accounting for the expected contribution 
of Poisson countin g unc e rtai nties. This is similar to the 
method that BeU et all (|20"07, ) applied to measure the 
dumpiness of the Milky Way halo. Although the reader 



is referred to this paper for a detailed description, the 
gist of the method is to compare the actual deviations 
between the data and the model with the expected de- 
viations from Poisson statistics. For a binned (but not 
smoothed) map of the considered system with N pixels 
of interest we therefore compute: 



V total/ 

TV 



N 



N 



Mi 



N 



4=1 



1 N 

4=1 



where Di is the number counts in the data for pixel i. Mi 
the expected number of stars in this pixel from the best 
model determined in Table 1 and Pi that of a Poisson 
realization of this best model. To avoid being sensitive 
to dumpiness in the background, induced e.g. by the 
clustering of misidentified background galaxies, we only 
consider the N pixels for which Mi is at least twice the 
number of st ars expected from the background alone. 
In contrast to iBell et all ()2007[ ) , who have a large num- 
ber of pixels, the value of N remains very low in the 
case of the maps of Milky Way satellites. To still get a 
good handle on the influence of Poisson noise, we rely 
on a Monte-Carlo analysis and generate numerous Pois- 
son realizations of the models to study the distribution 
of ((Tsc/total)^ values for each system. Moreover, since 
these values can become negative in the case of an ab- 
sence of distortion, we refrain from determining /total 
right away. In this case, (cTgc/total)^ < means that the 
scatter in the SDSS data is actually smaller than the typ- 
ical scatter exp ected froni shot-n oise for the best model. 

As noted bv iBell et all (|2007D . this method IS insen- 
sitive to the chosen pixel scale but will only be able to 
measure the influence of deviations from the best model 
that have a size at least bigger than the chosen pixel size. 
We therefore assume 4 different grids with pixel scales of 
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Fig. 4. — Smoothed maps of Bool from the SDSS data (left) and four simulations of Bool constructed assuming the best parameters 
of Table 1. The smoothing kernel is a two-dimensional Gaussian of dispersion 4'. The contours correspond to detections of 3cr, Scr and 
IOct over the background level, measured in an annulus between 0.75° and 0.9°. Although Bool has a distorted morphology, the purely 
spheroidal simulated models show as much distortion. 
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Fig. 5. — Distributions of (cTsc/total)'^ for the three dwarf galaxies Bool (left), CVnl (middle) and UMall (right). In each case, the 
results are shown for four grids with pixel sizes of 0.5' (left-dashed histogram), 1' (right-hashed histogram), 2' (white histogram) and 4' 
(grey histogram). Only the distributions obtained for CVnl with the 0.5' grid and for UMall for the 1.0' and 4.0' grids are different from 
at least at the Icr-level. The distributions obtained for Bool are typical of those obtained for the other dwarf galaxies. 



0.5', 1', 2' and 4' for which we determine the distribu- 
tion of (cTsc/total)^ values^^. The derived distributions 
are shown in Figure [5] for three examples: Bool where 
the pixels correspond to physical sizes of 10, 19, 38 and 
77 pc, CVnl (respectively 32, 63, 127 and 254 pc) and 
UMall (respectively 4, 9, 17 and 35 pc). The influence of 
the pixel grid is readily visible: the distributions derived 
from the small pixel grids have much more scatter, while 
larger pixel grids lead to almost no scatter and are closer 
to zero. In that case, there are not enough pixels for such 
sparse systems (and these examples are among the most 
populated of the new MW satellites). 

In most cases the distributions are similar to the ones 
of Bool shown on the left panel and imply the absence 
of any significant deviation from a scatter in the stellar 
distribution that is only produced by shot-noise. How- 
ever, CVnl and UMall are the two cases where Pois- 
son scatter alone may not account for the measured 
scatter in the distributions. The distribution from the 
0.5' = 32 pc grid of CVnl yields (crsc/total)^ = QM^oil, 
a mild la detection that might be induced by devia- 
tions from the best ML model in the data. There is 

It should be noted, however, that for close objects, the grid 
with the smallest pixels is not very reliable as most pixels are ex- 
pected to contain either or 1 star. This grid is only useful for the 
distant (D > 50 kpc) and populated objects in our sample. 



a slightly more significant deviation in UMall with the 
1.0' = 9pc and 4.0' = 35 pc grids yielding respectively 
(crsc/total)2 = 0.27ta\7 and (crsc/total)^ = 0.06 ± 0.03, 
or a contribution of intrinsic dumpiness to the scatter of 
52tll% and 24 ± 7% (the resuhs of the 2' grid, 22^^^%, 
although consistent with no distortion, have large enough 
uncertainties that they are also consistent with a signal 
of the order of that found for the 1' and 4' grids). 

It is reassuring that the two dwarfs that show the high- 
est evidence of dumpiness with this method are the two 
cases where distortions or substructure have been argued 
on other grounds. Indeed, a clump of young stars was un- 
covered in deep LBT observations of CVnl (jMartin et al.l 
HOOS) and, with a size of ~ 1', it should only be picked 
up by the grid with the smallest pixel size. In addition, 
UMall, the other galaxy that appears to contain a devi- 
ation from its best model is the likely progenitor of the 
Orphan Stream and is therefore expected to b e suffering 
from tidal disrupt ion (Fd lhaucr ct al. 2007; M artin et al.l 
120071 ahhough see lSales et al.ll2008f l. On the other hand, 
these two cases are only 1 — 2a deviations from null- 
detections and, given the 60 distributions generated for 
the 15 studied satellites, it cannot be ruled out that they 
are in fact chance detections that have nothing to do 
with real distortions in the morphology of the dwarfs. 
Although tentative, the detections in CVnl and UMall 
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therefore need to be confirmed from deeper, more popu- 
lated datasets. In all the other satellites, the scatter in 
the stellar density can be entirely accounted for by shot- 
noise in the SDSS data. Either these systems intrinsically 
contain no distortion or the SDSS is too shallow to allow 
detection of any such dumpiness. 

5. DISCUSSION 

We have re-derived the structural parameters and mag- 
nitudes for all the newly discovered Milky Way satellites 
(Table 1). Not only are these determined in a homo- 
geneous manner for all objects, but the ML algorithm 
produces robust results, even for extremely faint systems 
that only contain a few tens of stars in the SDSS. As men- 
tioned above, our results are overall in good agreement 
with previous measurements although there are some 
systematic differences, such as a higher ellipticity found 
with the ML technique compared to that measured from 
smoothed maps. Those cases are easily understood since 
sparsely populated smoothed maps tend to be dominated 
by the (round) smoothing kernel. 

The complete dataset can be summarized on the 
vs. My plane where dwarf galaxies and globular clus- 
ters were well separated prior to recent discoveries. Fig- 
ure [5] shows the Milky Way sat ellite system with globular 
clusters taken from the list of iMackev fc van den Bergtil 
(|2b05i) (with t he addition of the two extremely dim clus- 
ters found by iKoposov et al.l l2007bD represented as big 
dots, previously known dwarf spheroi dal galaxies as 
filled squares (the values are taken from lMateolll998D and 
recently discovered SDSS satellites as circled squares. All 
the new satellites are smoothly distributed from the re- 
gion inhabited by the brighter dwarf galaxies, down to 
ex tremely small and faint syst ems. As previously n oted 
bv lBelokurov et all (|2007f) and lGilmore et all (|2007f) . the 
clear size difference that existed before between globular 
clusters, with sizes smaller than ^ 20 pc and dwarf galax- 
ies > 100 pc becomes blurred at low luminosities. Many 
SDSS satellites are in fact smaller that 100 pc (BooII, 
CVnII, Com, SDSSJ1058-h2843, Segl and Will). There- 
fore, it might appear more difficult to assign these sys- 
tems to one population or the other. However, spec- 
troscopic surveys con ducted so far (jMartin et al.l 120071 : 
iSimon fc Gehal l2007f ) show that they are dark mat- 
ter dominated (except maybe for SDSSJ1058-I-2843 and 
Segl for which no spectroscopic data have yet been pub- 
hshed). 

Another obvious feature of Figure [5] is the higher el- 
lipticity of dwarf galaxies compared to globular clusters. 
Unfortunately, most faint G Cs (My > —6.0) have no 
ellipti city measurement in the IMackev fc van den Berglil 
(|2005f ) list. They also either lie away from the region 
covered by the SDSS or have incomplete photometry due 
to crowding at their center. As a result, we cannot de- 
termine their structural parameters using the ML algo- 
rithm. Still, in cases where crowding leads to incomplete 

These include Carina (Car), Dra, Fornax (For), Leo I, Leo II, 
Sculptor (Scu), Sextans (Sex) and Ursa Minor (UMi). The Sagit- 
tarius galaxy (Sgr) was removed from the sample as it is clearly 
interacting with the Milky Way and its properties are quite uncer- 
tain. As the only transitional galaxy of the sample, LeoT could 
arguably be removed from the following analyses, but that would 
not change our results given the mild (and uncertain) ellipticity of 
this dwarf. 



data in the SDSS catalog, the clusters do appear visually 
round in their outskirts. Pal 5 (My — —5.2) is one of 
the few GCs to which we can apply the ML algorithm. 
This object is of particular interest as it is known to be 
tidally disrupting in the gravity field of the Milky Way 
and shows beautiful stellar tidal tails (lOdenkirchen et a.1.1 
I2001bl: iRockosi et al.ll200l iGrillmair fc Dionatosll200ir 
Even in such an extreme case when one could expect the 
measured ellipticity to be inflated by the tails, applying 
the ML algorithm yields e = 0.11 ± 0.04 and the cluster 
remains very round. Moreover the 5 faint GCs in the list 
that have an ellipticity measurement all have e < 0.25 
so there does appear to be a significant difference be- 
tween the flattening of genuine GCs and that of newly 
discovered MW satellites. 

The left panel of FigurejTjshows the distribution of MW 
sa tellites in the e — My plane. As was already reported 
bv lvan den Berghl ()2008l ) using the previous, rounder el- 
lipticity estimates, GCs and dwarf galaxies populate dif- 
ferent regions of the plot (a Kolmogorov-Smirnov (KS) 
test actually yields a probability of 0% that both datasets 
follow the same ellipticity probability distribution) . Such 
a drastic difference is probably related to the formation 
mechanism of the objects in the two samples and/or to 
the presence of a dark matter halo around dwarf galax- 
ies whereas GCs are not believed to contain any dark 
matter. As for the dwarf galaxies, they appear to be- 
come more flattened as they become fainter. When di- 
viding the total sample of dwarf galaxies at a My — —7.5 
threshold^'^, more than a third of the faint systems (4 
out of 11) have ellipticities that arc higher than 0.5. On 
the other hand, this is the case for only a single bright 
satellite (UMi; or 2 if one accounts for Sgr). Further- 
more, three of these highly flattened faint systems (Her, 
UMal, UMall) are reasonably well populated systems in 
the SDSS {N.^. = 128, 70 and 296 stars respectively) and 
with only small uncertainties on e. 

To quantify this difference in the ellipticity, the right 
panel of Figure [7] presents the probability distribution 
functions of the faint (full line) and bright subsamples 
(dashed line). The arrows correspond to the mean of the 
two distributions, that are significantly different (a 4cr 
difference): e = 0.32 ± 0.02 and e = 0.47 ± 0.03, respec- 
tively (the uncertainties are determined from a Monte- 
Carlo resampling of the subsamples, redrawing the indi- 
vidual e values following the uncertainties of Table 1). 
A KS test also yields a very low probability (0.4%) that 
the two subsamples are derived from the same underly- 
ing distribution. The scatter of the PDFs is, however, 
hard to constrain as it results from two effects: from 
the larger uncertainties of faint systems, especially the 
rounder ones (e < 0.5), but also from the presence of 
objects whose uncertainties span a large range in ellip- 
ticities. In contrast, the majority of bright satellites are 
clustered in the range e — 0.25 — 0.35 (Car, Dra, For, 
Scu, Sex, LeoT). 

This somewhat arbitrary threshold was used as it divides the 
complete MW dwarf galaxy sample in two subsamples of similar 
sizes (10 and 11 members). This subsampling, that we do not 
claim to be physical in any way, also separates galaxies without 
any clear evide nce of complex ste l lar populations from those which 
do, hke LeoT (I rwin et aLll2007l : Ide Jong et"an I2008al) and CVnl 
f,Martin et al..i2008l ^! 
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rh (pc) 

Fig. 6. — Distribution of the MW satellites in the r^-My space. Colors code their ellipticity as defined on the wedge and blac k is used 
when no ellipticity measurement is available. Globular clusters, taken from Mackcy & van den Bergh (200D and IKoposov et al.l l|2007b), 
are shown as big dots and are rather small systems clustered on the left side of the plots; previously known, bright, dwarf galaxies are shown 
as filled squares (Mateo 1998), while new MW satellites appear as circled squares and become fainter and smaller than the other dwarfs, 
seemingly bridging the size gap that appears at higher luminosity between GCs and dwarf galaxies. Error bars represent the uncertainties 
taken from Tabl e 1 or the literature. The region of large and faint objects (bottom right) is incomplete due to surface brightness limits 
from the SDSS IKoposov et al.ll2007al) . 




M„ dP/dE 

Fig. 7. — Evolution of the ellipticity of MW satellites as a func- 
tion of their absolute magnitude (left panel). Globular clusters 
are this time shown as small hollow circles to avoid cluttering the 
plot, while the symbols for dwarf galaxies are the same as in Fig- 
ure m Error bars c orrespond to the uncertainties listed in Table 1 
or in IMateol 111998) for the brighter galaxies. Dividing the com- 
plete sample at My = — 7.5 (dotted line) leads to the probability 
distribution function (PDF) of the faint subsample, represented by 
a full line in the right panel, that extends to higher e values than 
that of bright galaxies (dashed line) . The mean of the PDF of faint 
(bright) systems is represented by the thick (thin) arrow. 



5.1. Why are the faintest satellites more flattened? 

Several possible explanations for these flat dwarf galax- 
ies have already been laid out in IColeman et al.l ()2007f ) 
for the case of Her, but are worth recalling and develop- 
ing in a broader perspective here. 

5.1.1. Could the faintest satellites be flattened by rotation? 

Current spectroscopic surveys, albeit usually targeting 
only a few tens of stars in each satellite, show no obvious 
sign of rotation. Even if the systems were rotationally 
supported, they would be approximately oblate. As we 
can assume we are seeing them at random viewing an- 
gles, their intrinsic flattening would be much more than 
< e >= 0.47 ± 0.03. In that case, they would have im- 
plausibly flat shapes and implausibly low intrinsic veloc- 
ity dispersions, given their already low measured velocity 
dispe rsions in the rangc4— 12kms~^ |K lc yna_et al, 2()05|; 
Muiioz et all l2006i : iMartin et all l2007l : ISimon fc Gehal 

5.1.2. Could the stars in the faintest satellites trace the 
shape of more elongated dark matter halos? 

It has been shown that dissipative gas infall dur- 
ing galaxy formation can induce changes in the shape 
of dark matter (DM) halos and make the m rounder 
than in the absence of any barvons fsee e.g. iDubinskil 
11994 iKazantzidis et al.ll2004l : iMaccio et al.ll2007t ). The 
lower baryon/star content of the recently discovered 
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MW satellites (iKlevna et all 120051: iMuiioz et all 120061 : 
iMartin et aH l2007t iSimon fc Gehal 120071 ) would permit 
the dark matter halo to remain triaxial, often nearly pro- 
late, for these systems. Assuming stars trace the DM 
potential, they would then appear more elongated in the 
f aintest systems . 

iKuhlen et all ()2007[ ) have studied the shape of the DM 
subhalos in the Via Lactea simulation of a Milky Way 
size ACDM halo. They follow subhalos down to a mass 
of 10® Mq, not much higher than those of the faint MW 
dwarf galaxies studied here, and find them to have a more 
elliptical shape in their central regions. This could pro- 
vide an alternate explanation for the shape of the faintest 
satellites, which are also the smal lest (because the large r 
ones cannot be found in the SDSS: IKoi30SOv et al.ll2007al ): 
stars probe a smaller region of the subhalo, and therefore 
trace out a less spherical potential. This latter argument 
cannot, however, explain the shape of the three satellites 
with the highest e (Her, UMal and UMall) as all three 
are large satellites with > 140 pc. 

The main issue with this scenario is the assumption of 
stars distributing themselves in the DM subhalo so as to 
follow the shape of underlying mass potential for which 
no cogent mechanism exists. 

5.1.3. Could the elongation be the sign of tidal distortion? 

The best example of tidally induced elongation in 
the Local Group is given by the Sagittarius dwarf that 
is clearly being disrupted by the Milky Way as it 
leaves behi nd stellar stream s wrapped ar ound the Galaxy 
(jlbata et a l. 2001; Maicw ski et al.ll2003f ). This resuhs in 
an ellipticity of ~ 0. 65 for the (already dis rupting) main 
body of the dwarf (jMajewski et al.l[2003f ). close to the 
values that we have measured for Her, UMall and, to a 
lesser extent UMal. The absence of any direct evidence of 
stellar streams originating from the faint dwarf galaxies 
could easily be explained by their faintness that implies 
that, if they exist, such features should be very sparsely 
populated. UMall has also been proposed to be the pro- 
genitor of the Orphan Stream, which would lend support 
to this tidal distor tion scenario, a t least for this satellite. 

As explained in Coleman et al.l (|2007[ ). tidal distortion 
requires the satellite to pass close enough to the MW 
to be stripped. Simple arguments can be used to es- 
timate the required pericenter of the galaxy's orbit as 

^peri - '■gal ' ^Mw/CTgal with Tgal itS sizC SCale, CTgal itS 

velocity dispersion and ctmw — 150 km s^^ that of the 
Milky Way. For Her, UMal and UMaH, the inferred peri- 
center is of the order of ^ 10 kpc which makes for very 
eccentric orbits in the case of UMal (i?apo > 97 kpc) and 
Her (i?apo > 132 kpc), with an eccentricity of e > 0.90 
for both. In addition, this latter object is certainly 
not at apogalacticon since it shows a large receding 
motion from the MW {vr — 145 km s^^, corrected for 
the Solar motion). Very eccentric orbits are not un- 
heard of in dark matter simulations of galactic groups. 
For example. Gill et al.l (j2004l ) show that their surviv- 
ing satellites at z = have a mean orbital eccentricity of 
< e >= 1 — RpprJ Rnpn = 0.61, with a standard deviation 
of 0.19. IBensonI ()2005[ ) also finds that half of infalhng 
dark-matter substructures are on near parab olic orbits 
(e = 1.0 ± 0.1). Moreover, 'Sales_eLalJ (|2007t ) identified 
in their simulations of galaxies that about a third of 
their satellite dark subhalos have 'unorthodox' (includ- 



ing high eccentricity) orbits from three body interactions 
in the galactic halo. In this scenario, the satellites most 
severely stripped by tides would be both the least lumi- 
nous and th e most elongated. 

However, iPenarrubia et al] ()20 08b) have recently ar- 
gued that the observed size and central radial velocity 
dispersions of the new SDSS satellites are generally not 
low enough to be disrupted versions of the Fornax, Draco 
or Sagittarius dwarfs. So this explanatory scenario also 
appears problematic. 

Overall, none of these three scenarios is entirely satis- 
factory but we believe that the latest, elongation induced 
by tidal interaction with the Milky Way, may be the least 
problematic. Indeed, there is a known stellar stream that 
could be produced by the disruption of UMall, the Or- 
phan Stream, and a few MW satellites with very eccentric 
orbits is far from impossible. M oreover, the importance 
of the lPeharrubia et al.l ()2008bl ) caveat depends, for each 
satellite, on the amount of stellar loss, which is all but 
constrained. The new satellites could therefore corre- 
spond to a murky sample of galaxies at different stages 
of their disruption, strongly hampering any analysis of 
the population as a whole. 

6. CONCLUSION 

We have presented a rigorous and homogeneous (re- 
) analysis of the structure of the stellar bodies of the re- 
cently discovered Milky Way satellites. The application 
of a Maximum Likelihood algorithm to SDSS data yields 
robust estimates of the centroid, position angle, elliptic- 
ity, exponential half-light radius and number of mem- 
ber stars (within the SDSS) for the 15 faintest candidate 
dwarf galaxies around the Galaxy. With the assumption 
of stellar population models, derived through the analy- 
sis of the same dataset, we have also studied the impact 
of 'CMD shot-noise' on the observed luminosities and 
stellar masses of these satellites. The magnitude of the 
faintest systems can be hampered by a large scatter (of 
up to ^ 1 mag) that is due solely to the individual evolu- 
tion of their very few stars. Stellar mass estimates suffer 
less from this sparseness but, obviously, depend strongly 
on the assumed IMF. Finally, we have shown that except 
for tentative evidence of deviations from the best model 
in the CVnl and UMall dwarf galaxies, the apparently 
distorted morphology of the satellites in the SDSS can 
be entirely accounted for by shot-noise. This does not 
mean that some distortion is not present, but only that 
deeper datasets are required: one should not rely on a 
visual inspection of stellar density maps as the structure 
in these is strongly driven by the low number of stars 
that are observed in each satellite. 

In our discussion, we focussed on the remarkable re- 
sult that the faintest satellites are, on average, highly 
flattened. We put forth three lines of explanations, 
but none of these three scenarios is entirely satisfac- 
tory. However, we believe that tidally induced elonga- 
tion may be the least problematic and is certainly the 
easiest to (dis)prove from observed distortions or tidal 
streams from the stellar maps of these dwarf galaxies. 
This however requires deeper follow-up data than the 
SDSS since our search of such distortions yields mainly 
null-detections when accounting for Poisson noise. Alter- 
natively, a final answer will strongly benefit from deeper 
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systematic su rveys such as the PanSTARRS Stt survey 
(|Kaiseill2006l ) that should significantly expand the SDSS 
harvest of faint satellites of the Local Group. 
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APPENDIX 

TWO-DIMENSIONAL RELATIVE LIKELIHOOD DISTRIBUTIONS 

The following contour plots show, for each satellite, the evolution of the relative likelihood for all the combinations 
of two parameters that are fit by the ML algorithm. In each panel, outward from the best model that is represented 
by the dot, contours represent drops of 50%, 90% and 99% of the likehhood with respect to the best model. 
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Fig. 8. — Evolution of the relative likelihood for all the combinations of two parameters that are fit by the ML algorithm for Bool, BooII, 
CVnl and CVnII. 
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Fig. 9. — Same as Figure[8]for Com, Her, LeoIV and LeoT. 
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Fig. 10.— Same as Figure|8]for Segl, UMal, UMall and Will. 
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Fig. 11.— Same as Figure|8]for SDSSJ1058+2843, Object Y and Dra. 
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